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Abstract 

We study the conformal boundary conditions of the dilute 0{n) model in two dimensions. 
A pair of mutually dual solutions to the boundary Yang-Baxter equations are found. They 
describe anisotropic special transitions, and can be interpreted in terms of symmetry breaking 
interactions in the 0{n) model. We identify the corresponding boundary condition changing 
operators, Virasoro characters, and conformally invariant partition functions. We compute 
the entropies of the conformal boundary states, and organize the flows between the various 
boundary critical points in a consistent phase diagram. The operators responsible for the 
various flows are identifled. Finally, we discuss the relation to open boundary conditions in the 
0{n) model, and present new crossing probabilities for Ising domain walls. 

1 Introduction 

Loop models in statistical mechanics have been studied for a very long time, both for their physical 
and mathematical properties. They have enjoyed considerable interest recently as fundamental 
objects in the SLE approach to critical phenomena, and in the search for systems with quasiparticles 
obeying non-abelian statistics [UOISj. 

Although a lot is known about these models, new, very basic features keep being discovered. An 
example of this concerns conformal boundary conditions in the two dimensional (classical) case. 

Conformal boundary conditions (CBCs) are an important tool in the understanding of conformal 
field theories (CFTs) and in many applications. They are systematically classified for minimal 
models, and usually complicated to implement physically, except for the simplest ones: "free" and 
"fixed", in the proper variable. Surprisingly, it turns out that for dense loop models (see the 
definition below), a continuous set of CBCs is obtained simply by deciding to give loops that touch 
the boundary a fugacity ni different from the fugacity n for loops in the bulk. These CBCs — 
now dubbed "JS boundary conditions" — have been studied in a series of papers using lattice and 
algebraic techniques [H [U [H [7] • They have also been tackled in terms of matrix models f8l|9], and 
found applications in the study of logarithmic CFTs \W\ lllj. 

There are basically two kinds of loop models: the loops are either dilute or dense (for recent 
reviews, see [l2l SS]). In the former case, the loops occupy only a critical fraction of the available 
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lattice sites, while in the latter case, the loops cover more than this critical fraction. Bulk properties 
then are the same, independently of this fraction, and characterized by what is called the dense 
universality class (often, this universality class is most conveniently studied on the square lattice by 
forcing the loops to cover the whole lattice, giving rise to a fully packed loop (FPL) model. One has 
to be careful however that properties in the FPL case are not always universal) . If n is the fugacity 
of the loops, the dilute theory is believed to be in the same universality class as the 0{n) spin 
model (once continued analytically to real values of n), and conformal for n g [—2,2]. The dense 
theory on the other hand does not describe the low temperature phase of the generic 0{n) model, 
as self-intersections are absent, which now play a crucial role [14j. The classification of conformal 
boundary conditions for 0{n) models in their Goldstone phase is still an open problem. 

The dilute model admits one more relevant parameter than the dense model: the fugacity of 
the basic monomers making up the loops. Only when this parameter x is adjusted to a critical 
value Xc does the model become conformally invariant in the bulk. This freedom translates into a 
more complicated boundary behaviour: not only can one adjust the weight of loops touching the 
boundary, but one can (and must) also adjust the weight of boundary monomers. Indeed, in the 
dilute model a generic loop touches the surface with probability zero, so the surface fugacity has to 
be critically enhanced to allow a finite fraction of boundary sites to be occupied by loops. This is not 
necessary in the dense model since the boundary is always covered with loops in the continuum limit. 
This has also a nice interpretation in the language of symmetry breaking boundary interactions in 
the 0{n) model [15], that we discuss shortly below, and is the key to physical appHcations to appear 
elsewhere [16J. Once this feature is under control, the basic aspects of CBCs in the dilute case are 
formally similar to those in the dense case, after a proper redefinition of the parameters. Carrying 
this out in details gives us control of key combinatorial quantities and crossing probabilities that 
were not known up to now. 

The paper is organized as follows. In section 2 we introduce the problem in more details and 
discuss the basic features of the phase diagram. In section 3 we find a new solution to the boundary 
Yang-Baxter equation, which generalizes the solution of [17] in the special case to what we call the 
"anisotropic special case". This constitutes the new CBCs for the dilute loop model. In section 4 we 
discuss the basic features of these new CBCs, and determine the associated critical exponents and 
annulus partition functions. In section 5, we discuss some features of the boundary phase diagram, 
in particular the existence of RG fiows and their relations to the boundary entropy. In section 6, 
we discuss yet another type of CBC, where, instead of affecting the weights of loops touching the 
boundary, we introduce a boundary magnetic field, and thus "open" the loops on the boundary. In 
the dense case, it turned out that these open boundary conditions could be reformulated exactly 
as a particular case of the JS boundary conditions. The situation is a bit more complicated here, 
however similar. We finally apply our results to the determination of crossing probabilities of Ising 
clusters on an annulus. In section 7 we give a more detailed summary of our results, and we discuss 
some further generalizations and open directions. 

2 Special transitions in the 0{n) model 

We start by giving a small review of known results about the 0{n) model and its (bulk and surface) 
critical behaviour. We do not try to be exhaustive, and refer the interested reader to the vast 
literature on the subject [111 [13l [ill [19l [20l [HI [22] for further details. Nevertheless, we want to set 
up a self-contained physical framework in which further discussions will appear most natural. 



2 



2.1 The 0{n) model and its loop expansion on the honeycomb lattice 

The 0{n) model is usually presented as follows. It can be defined on an arbitrary lattice in D 
dimensions, although we will soon restrict ourselves to D = 2 and to the honeycomb lattice. It is 
a spin model, with n-components spins S'f Hving on the sites i of the lattice, with fj. G {1, . . . ,n}. 
The spins are subject to the constraint Si. Si — ^i^i — which means they actually live on 

a (n — l)-dimensional sphere. The partition function of the model is taker0 as 




where the (ij) are the lattice edges in the bulk, x is a parameter interpreted as an inverse tem- 
perature. In particular, note that the C(n = 1) model is the Ising model with x — tanhiiT if the 
partition function of the Ising model is written Zising — ^ Y\. e^^'^' ■ The trace Tr {.} is pro- 

conf. {S} <ij> 

portional to J dSiS{Si.Si — n). It can be normalized such that it has the properties Tr {1} = 1, 

i 

TrISfS]'} andTr{Sf} = Tr{(Sn'} = Tr{(Sf)'} = ... = 0. 

With these relations, the partition function H]) allows a "loop" expansion. This loop formulation 
is well-known, and is obtained by a high-temperature (small x) expansion. When we take the trace, 
the terms with odd powers of S^ cancel. The terms of the form S^^ S^^^ S^^ S^^ . . . S^^ S^^^ subsist, 
where {1112), (*2*3), etc. and {ikiij are links of the lattice. Each such term can be represented 
as a loop on the lattice, and gets a weight n once the trace is taken. Other terms with an even 
power of 5^ , such as (Sf ) can be nonzero. Such terms can be avoided if we restrict ourselves to 
trivalent lattices (with at most three edges connected to each site). This is the reason why the 
two-dimensional 0(n) model is usually defined on the hexagonal lattice. Of course, one could argue 
on universality grounds that in the end the critical behaviour of the model is lattice independent. In 
any case, for the honeycomb lattice the loop formulation is exact, and we end up with the celebrated 
loop partition function 

conf. 

where the sum is taken over all the configurations of non-intersecting loops that one can draw on 
the honeycomb lattice. X is the number of monomers for each configuration, while N is the number 
of loops. Note that the loop model with the partition function ^ is more general than the spin 
model (HJ), because in the former case n is not restricted to be a positive integer. However, both 
interpretations of the 0{n) model, as a magnetism model or as a geometric loop model, are useful 
to get some intuition about its rich critical behaviour. 

^Of course, a more natural 0{n) model would be Z = Tr |n<ij> ^^P (^•^i 'Sj) }, and this is indeed the generic 
0{n) model which has been studied for a very long time and in higher dimensions. When D = 2 however, one usually 
prefers to consider the model ([TJ as it is much simpler and yet it is believed to belong to the same universality class 
up to the critical point. In the Goldstone phase however, the two models are believed to have different critical 
behaviour generically (the Ising case n = 1 is an noteworthy exception), as discussed briefly in the introduction (see 
also reference [f4)V 
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As it is a ferromagnetic spin model, we expect the 0{n) model to possess a disordered phase 
when X is small, and an ordered phase when x is big. In two dimensions, the transition between 
these two regimes is known to be of first order when n > 2. Therefore we restrict to — 2 < n < 2, 
where the model exhibits a second-order phase transition. In this case, the following phase diagram 
holds 

No loops Dilute loops Dense loops Fully packed loops 
• « • W • « • 

X = X = Xc X ^ Xo X ^ oo 



Massive theory Dense loops phase (massless) 

Disordered spins Ordered spins 

In terms of loops, this can be understood as follows. For a; = there is no loop at all. For small 
a: > 0, some small loops appear and as x increases, the loops are more and more numerous and can 
become longer and longer. At some point {x = Xc) the mean length of the loops going through a 
randomly chosen point diverges. If one increases the fugacity x a little more, the loops proliferate 
and they finally cover each edge of the lattice with finite probability (even when the lattice is 
infinite), so they fill the whole space in the continuum limit. Note that in the discrete setting 
however, even for x > Xc, the loops do not fill all the sites of the lattice. But when x = oo, that is 
in the fully packed loop (FPL) phase, all the sites of the lattice are covered by a piece of loop. This 
is a rather strong constraint on the loop configurations, and in particular the universality class of 
the FPL phase is not lattice independent. We will not be concerned by this theory in what follows, 
however it should be noted that it corresponds (on the honeycomb lattice) to a non-trivial theory, 
different from the one of the dense phase. The dilute phase {x = Xc), the dense phase {xc < x < oo) 
and the FPL phase {x — oo) are described by conformal field theories with central charges 



Cdiiuto = 1 - 6 g = 1 + - (3a) 

g TT 

cdcnso = 1 - 6 ^-^ / 5' ^ 1 _ I (3b) 

9 71" 

CFPL = 1 + Cdonse (3c) 

where n = 2 cos 7 with 7 e [0, n). Note the consistency with Zamolodchikov's c-theorem ^2^, which 
states that the central charge always decreases along the RG 

The exact localization of the points xq and Xc for the honeycomb lattice is the point of Nienhuis' 
conjecture [19] 

v/2 + 72^ V2 - 

which we will discuss in greater detail in section [3l 



^This theorem is true for unitary theories, which is not the case in general for the 0{n) model. Therefore the 
application of this theorem here is not justified and certainly not rigorous, but it still provides some insight in the 
behaviour of the model. 
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At this point one might be worried by the special case n — I, which looks a bit different from 
the traditional results about the Ising model. However, it is not difficult to see that everything 
fits together here if we recall that x is related to the usual Ising coupling by a; = tanhif. Then 
X = xq = 1 corresponds to K = +00, and x > xq gives a complex couphng K. Hence, the zero 
temperature point of the Ising model is x = xq, and the part x > xq is irrelevant when we deal 
with the model with K real. Also, the Ising model is known to be dual to an Ising model on the 
dual lattice, so the theories aX K = and K = +00 should be dual to each other. Again there is 
no contradiction here, since Cdonso = 0, as expected if it has to be related to a massive theory. 

2.2 Surface critical behaviour : ordinary, extraordinary and special tran- 
sitions 

So far, we have focused on the bulk critical behaviour of the 0{n) model. Now let us turn to its 
behaviour near a boundary. We consider the partition function 




which allows the possibility to enhance the coupling y between two boundary spins above its bulk 
value X. To begin with, let us drop this possibility and assume that the spin couphng is the 
same for the spins living on the boundary and those in the bulk: y — x. When one goes from 
X < Xc to X > Xc, the bulk spins order. However, the boundary spins have less close neighbours 
and therefore the spontaneous magnetization vanishes near the boundary, on a typical length of the 
order of the bulk correlation length. This transition to ordered spins in the bulk, with a spontaneous 
magnetization in the bulk which drops to zero near the boundary, is called the ordinary transition. 
Now, consider x < Xc'- the bulk spins are disordered. If the boundary coupling y is sufficiently 
enhanced above x, then the boundary spins order, despite the fact that the bulk is disordered. This 
transition from disordered to ordered spins in a one-dimensional layer along the boundary of width 
of order of the bulk correlation length, while the bulk remains disordered, is called the surface 
transition. It is expected to belong to the universality class of the one-dimensional 0(n) modeH. 
If we start from a configuration with the surface spins ordered and the bulk disordered (x < Xc), 
and go to X > Xc so that the bulk spins order, the transition is called extraordinary transition. 
If we follow the surface transition line y{x) when x ^ Xc, the boundary and the bulk correlation 
lengths both diverge, and this point is known as the special transition. At this point, the coupling 
y is critically enhanced, in the sense that bulk and surface order simultaneously. In other words, 
y precisely compensates for the lack of nearest neighbours at the boundary. This bulk/boundary 
critical behaviour is summarized in figure [H 

^Of course there can be such a boundary transition at finite coupling K only if the lower critical dimension Dlc 
of the 0{n) model is such that -D_lc < 1. In particular, this would imply n < 1. However, as in the case of the Ising 
model above, we will actually never work with the real coupling K or its boundary equivalent Ki,, but rather with 
X = tanhK and y = tanhKi,. Thus nothing should prevent us from considering situations where y > I, even if it 
does not correspond to a real coupling Kc. In the loop model, the parameters x and y are the natural parameters, 
and the surface transition exists in this model at finite y. Thus in this paper we will always avoid this discussion 
about the lower critical dimension, and we will not restrict ourselves to n < 1 (as the authors of j22j did), because 
even if the surface and extraordinary transition do not appear in the 0{n) spin model at finite Ki,, they exist as 
geometric transitions in the loop formulation of the model. We keep discussing these transitions in the spin language 
because it still gives some intuition about the critical behaviour. 
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Figure 1: Phase diagram of the 0{n) model with bulk spin coupling x and boundary spin coupling 
y. In this paper we will focus on the line x = Xc, which corresponds to the dilute phase for the 
loops. 

The loop interpretation of the special transition goes as follows. The partition function ^ 
becomes 

z=J2 ^''y'^^'' (6) 

conf. 

where N is the number of loops, X is the number of bulk monomers, and Y is the number of 
monomers at the boundary (for each configuration). When we focus on the line a; = Xc, as we will do 
in the rest of this article, we have a model of dilute loops with fugacity n which get different energies 
per unit length when they come to touch the boundary. Since we now work with the following one- 
dimensional phase diagram, there is no point in talking about ordinary or extraordinary transitions, 
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and we will adopt the terminology of ordinary and extraordinary boundary conditions (b.c) instead. 

Special b.c 

• • *-* • *^ • 

y = y y^ys y ^ oo 



Ordinary b.c Extraordinary b.c 

When y is small, the loops avoid the boundary, and behave as if they lived on a lattice with one less 
row on which we put free boundary conditions. Hence the theories with y = or y = x = Xc are 
obviously the sam^, and correspond to free boundary conditions for the loops. On the contrary, if 
y is big enough, then there is one polymer which is adsorbed along the boundary, and it prevents 
the other loops from touching the surface. These loops are then not affected by the surface, and 
they behave as in the case of ordinary b.c. Hence the RG fixed points describing the ordinary and 
extraordinary b.c correspond to the same theory, except for the fact that there is one adsorbed 
polymer on the surface in the extraordinary case. This is illustrated in figure [2l 

The transition between the free regime and the adsorbed one is the special transition. The special 
point corresponds to a non-trivial scaling limit, where the loops can touch the boundary without 
being completely glued on it. To our knowledge, the exact localization of the special point on the 
honeycomb lattice has first been conjectured by Batchelor & Yung [17] 

2/5 = (2 - ny'/^ (7) 
and again, we will come back to this result in section [31 




Figure 2: The ordinary (a) and extraordinary (b) boundary conditions for the dilute loop model. 
The ordinary b.c corresponds to a free b.c for the loops. With extraordinary b.c, one loop is 
adsorbed on the surface, preventing the remaining loops to touch it. The latter then behave exactly 
as if there was a free b.c. 



^That is the reason why we do not indicate a RG flow between the points y = and y = Xc. 
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We would like to conclude this section with a remark about the dense phase of the 0{n) loop 
model. For x > Xc one could also consider a model where y is enhanced, and wonder whether there 
is a special transition in the dense loop model or not. The answer is no. The reason for this is that 
in dense loop models, the boundary is already covered with loops, so one cannot change the critical 
behaviour by encouraging the loops to touch the boundary. We expect that as soon as y < cjo, the 
coupling renormalizes towards y = x. However, in the present model on the honeycomb lattice, the 
situation y = oo still corresponds to a loop adsorbed on the boundary, exactly as in the dilute phase 
{x = Xc). This means that, when we diagonaHze the transfer matrix of the model, one half-loop is 
adsorbed on the boundary while the other half-loop moves freely in the system. Then the effective 
central charge we compute should be Ceff = Cdcnsc — 24/ii 2 rather than Cdense (^1,2 being the one-arm 
exponent in the dense phase) . A numerical check of these observations is presented in figure [3l It 
is obtained by computing the free energy per site Jl for successive L = 8, 9, . . . , 15, and relating 
the finite-size corrections of /l to the effective central charge by the well-known relation 



up to order O (-jJj). More generally, the behavior discussed in [4],[5] was technically obtained only 
for loops covering the whole square lattice, that is a maximally dense situation. It turns out that, if 
one relaxes this constraint to the more general dense case, where loops cover more than the critical 
fraction of available lattice sites (that is, Xc < x < 00), the same behavior is obtained: the same 
formulas for conformal weights apply, independently of the fugacity of surface monomers, and of 
bulk monomers, that is, independently of x and y. 

2.3 Anisotropic special transition 

Following Diehl & Eisenriegler [15], we introduce an anisotropic generalization of ([5]). The idea is 
to break the 0{n) symmetry of the interaction at the boundary down to 0(ni) x 0{n — ni), which 
allows the coupHng to be different for the ni first components of the spin and for the n — ni other 
ones. This can be written 



which of course makes sense only when n and ni are both integers. The loop framework, however, 
provides an analytic continuation to non-integer numbers of components ni, exactly as it did for 
the bulk parameter n. 

The boundary critical behaviour of this model was studied perturbatively by Diehl & Eisenriegler 
|15) . They used scaling arguments, mean-field theory and e-expansion below the upper critical 
dimension 13 = 4 to deduce the phase diagram and the crossover exponents appearing in their 
model. Let us sketch some of there results as follows. The case when y^ = yn-m is called 
isotropic and was discussed in section 12.21 when the bulk crosses the line x — Xc and goes from 
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Figure 3: Influence of the surface monomer fugacity y in the dilute (a) and dense phase (b) of the 
0{n) loop model. Here we have n — -\/2, x — Xc ^ 0.60 for (a), a; = a;o « 0.90 for (b). In (a) the 
CFT has central charge 7/10. We see that there are three regimes: y < ys where CcS — c = 0.7, 
y — ys with Ccff — c — 24/ii_2 — —1.7 and y > ys with CcS — c — 24ft,2,i — —9.8 (one-arm exponent 
in the dilute phase). In (b) the central charge is 1/2. Although the numerical precision is not good, 
our conclusion is that CcS = c — 0.5 for every flnite y (indeed, when the size N increases, Coff(A^) 
goes to c = 0.5). For y = oo it is Ccs = c — 24/ii 2 = — 1 (one-arm exponent in the dense phase). 

its disordered to its ordered phase, the boundary behaviour depends on the boundary coupling y 
and three different transitions are possible : ordinary, extraordinary or special. In the anisotropic 
case, things are quite similar. Let us start with j/„j and j/n-m both small. Clearly, this cannot be 
very different from the isotropic case when y is small : because of the lack of close neighbours, the 
(anisotropic) coupling between the boundary spins is weaker than in the bulk and the spontaneous 
magnetization cancels at the boundary. This means that in the space of parameters {ynnUn-m), 
the RG fixed point of the ordinary transition is attractive. 



In the opposite case, when j/„j and y„-ni are both large, the boundary spins are ordered even 
when the bulk is disordered, and there is a spontaneous magnetization at the boundary. This has 
already been discussed in the isotropic case, and corresponds to the extraordinary transition when 
X Xc- When ym > yn-m and j/„i, j/n-m 00, the boundary spins all point in the same 
direction S, where S''^ = for /3 G {ni + 1, . . . , n}. The dual case y„_ 
ground state with all the boundary spins having their components {1, 



> yni would lead to a 
, ni} all equal to zero. 
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Vn—ni 



Figure 4: Generic phase diagram for < ni < n in the rotated {yni^yn-m) plane. 

Crossing the isotropic half-hne ym = Vn-m > ys corresponds to switching from one ground state 
to another, and the transition is of the first order. 

Wiiat happens wiicn we go from the regime of tiic ordinary transition to the (anisotropic) ex- 
traordinary one? Let us fix yn-m to some value (smaller than ys). When y„^ is sufficiently small 
the theory fiows towards the attractive point corresponding to the ordinary transition, and when 
yni is sufficiently large it goes to the (anisotropic) extraordinary transition with the last n — ni 
components of the spins all equal to zero. The transition between the two regimes happens when 
the m first components of the boundary spins order. This transition is cafied the anisotropic special 
transition. For fixed values of n and ni, we can consider also the dual case y^ < yn-ni which 
leads to a different anisotropic transition. Then, in the space of parameters {yni^yn-m) there are 
three regions : two of them correspond to the (anisotropic) extraordinary transitions, and one to 
the ordinary transition. The three of them are separated by : the isotropic line between the two 
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extraordinary transitions (and this is a first-order transition) , and the two anisotropic special tran- 
sitions between the ordinary and extraordinary transitions. These three fines all meet at the special 
transition point Sp (figured]). 

RG fiows can be organized as fofiows. As argued above, there is a fixed point for the ordinary 
transition, which we cafi Ord, and this point is stable. There are two extraordinary points at 
infinity, which are stable. On each of the two special transition lines, there is one RG fixed point 
ASni or ASn-m ■ We will argue in the following parts of this article that these points are stable for 
a perturbation along the line, but unstable in the other direction. Note that, all along the line the 
theory is attracted by this point and hence is in its universality class. The special transition point 
Sp is unstable in both directions. This is all summarized in figure HI 




Figure 5: Two different kinds of loops can touch the boundary. They are marked with a circle 
(blob) or with a square. Blobbed loop (with circles) have a fugacity ni, and unblobbed ones (with 
a square) a fugacity n — ni. In the bulk, the loops still have a fugacity n. 



To discuss the anisotropic transitions in the the loop language, we have to introduce a sim- 
ple object which we use to distinguish some loops when they touch the boundary. This ob- 
ject is the 6/06 (see figured]). It is used to mark the loops coming from the trace of the terms 
^ii ■ ■ ■ ^ii ((*i*2), (*2*3), stc. and (ifcZi) are links of the lattice) restricted to the 

components a G {1, . . . , ni} in ([9]). These blobbed loops have a fugacity ni instead of n. Of course, 
there is an equivalent object for the n — m remaining components : these loops will be called the 
unblobbed ones, and have a fugacity n — m. Each time a blobbed loop touches the boundary it takes 
an additional blob. This does not change its fugacity, but the weight of the monomers touching the 
boundary is different if they are blobbed or unblobbed. To simplify a bit the formulae to come in 
what follows, let us introduce the weights 

^. ^ [^^y (lob) 
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which can be viewed as the fugacity of the circles and squares (figure [5]) when all the monomers get 
a weight x, without distinction for the boundary ones. The partition function ([9]) becomes 

Z=Y,x'' u;^* y^"^' <• (n - n,f' . (11) 

conf. 

In this formula, N is the number of bulk loops, N, is the number of blobbed loops and the 
number of unblobbed loops. X is the total number of monomers (boundary and bulk ones), W, 
is the number of blobs and W, is the number of "unblobs" (squares) . When we work within this 
loop picture, the phase diagram shown in figure H] must be interpreted as follows. Let us again fix 
a; = Xc, so we are dealing with a critical dilute loop model in the bulk, and we adopt the terminology 
of ordinary /extraordinary b.c instead of transitions. For n and ni fixed, there are three different 
regimes: the ordinary b.c when the loops avoid the boundary, and two (anisotropic) extraordinary 
b.c when there is one polymer (blobbed when w, > Wm, unblobbed when > w,) adsorbed all 
along the boundary (see figure [2]). The transition between the two extraordinary b.c is of the first 
order. The transition between the ordinary and one of the extraordinary b.c is the anisotropic 
special transition. Again, it corresponds to a non-trivial scaling limit, where there are two kinds of 
loops (blobbed and unblobbed) which touch the boundary without being adsorbed on it. We expect 
the geometric characteristics to be different for the blobbed and unblobbed loops. The universality 
class of this transition is completely characterized by the RG fixed point which sits on the transition 
line, and is stable under perturbations along the line. This point was named AS^ (or ASn-m) in 
the context of the spin model, let us rename it AS, (or ^45*,) for convenience when we deal with 
the loops. The two anisotropic transition lines and the first order transition meet at the (isotropic) 
special point Sp. We already discussed that point in the context of the dilute loop model in section 
12.21 In section [31 we will arrive to the following conjecture for the exact position of the points AS, 
and AS, for the model on the honeycomb lattice 



(AS,) { ' , ^j^^ ^ (12a) 

= 1 + iV2~?i 
l + iV2" 

In the next sections, we will be interested in the anisotropic special transition in a more quantitative 
way. 
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3 Integrable and critical points in 2D : from bulk to boundary 

Statistical models in two dimensions can be exactly solved when they are related to solutions of the 
Yang-Baxter equation. This equation implies the existence of a continuous family of commuting 
periodic transfer matrices (or Hamiltonians) which can be used to compute many quantities - such 
as the low-energy spectrum - using Bethe-ansatz techniques. In the case when there are boundaries, 
the equivalent of the Yang-Baxter equation is the Sklyanin refiection equation |25], sometimes called 
boundary Yang-Baxter equation. With a i?-matrix solution of the Yang-Baxter equation in the bulk 
and a if-matrix solution of Sklyanin's equation at the boundary, one can build a continuous family 
of transfer matrix with boundaries ^25] . 
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Although there is no general result relating integrability and criticality of a model, they can some- 
times coincide. Integrable points often play some particular role in the phase diagram of a model, 
and in general the information they provide is a key point to understand the critical behaviour. 
In the particular case of loop models, some signs of a deeper relation between integrability and 
criticality have appeared very recently in the literature [261 123 ■ The relation between integrability 
in the sense of Yang-Baxter and the lattice holomorphicity of certain discrete observables, which 
would lead to critical models in the continuum limit, is a fascinating subject. We hope that future 
work will develop this approach and provide a straightforward and rigorous way to relate lattice 
parameters to the characteristics of the objects showing up in the continuum limit [28]. But for now, 
let us use the more traditional way of studying a lattice model: we look for (bulk and boundary) 
integrable solutions, and conjecture that they correspond to critical points of the loop model. 

3.1 The loop 0{n) model on the square lattice 

As explained in part [21 the 0{n) spin model can be reformulated as a dilute loop model on the 
honeycomb lattice. To study the integrability of the model, it is useful to relax the constraint that 
it lives on the honeycomb lattice, and to define a dilute loop model on the square lattice [29]. The 
model contains 9 different plaquettes in the bulk. Two neighboring faces share a common edge, 
which can be crossed or not by a single loop. A loop cannot be stopped on an edge {ie such a 
configuration is given a weight zero). The total Boltzmann weight of a configuration is given by 
the product of the weights of the plaquettes and a factor n^, with N the number of loops. 




UJl 0J2 IJJ3 1^6 



Vertices related by a horizontal or vertical refiection are given the same weight, so there are only 
6 independent weights Wi, . . . , wg. An i?-matrix solution to the Yang-Baxter equation is given by 
the sum of the above 9 diagrams with weights [29j 

uji = sin2<i>sin3$ + sinusin(3$ — w) 

LU2 = sin 2<i> sin(3$ - u) 

= sin2<i>sinu , . 

UJ4 — sinusin(3$ — u) 

= sin(2<i> - u) sin(3$ - u) 

uje = — sinusin(<i> — u) 
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where u is the spectral parameter and $ is related to the weight n of a closed loop by the parame- 
terizatiorll 

n=-2cos4<I>. (14) 

It is convenient to draw the i?-matrix as 



R{u) = 




where it is understood that the marked angle stands for the orientation of the plaquettes. Note 
that, in particular, the weights (fT3)) imply that 







and that when the spectral parameter u goes to zero, the i?-matrix is proportional to the identity 



sin 2$ sin 3$ 





Other relations are satisfied by this i?-matrix, but we will not need them in what follows. We refer 
the interested reader to 1171 [29l for more information. 



We want to turn back to the honeycomb lattice, the present digression about the square lattice 
being only a trick to catch an integrable model (a solution to the Yang-Baxter equation). In the 
above model, the square lattice can be viewed as the honeycomb one for a specific choice of the 
spectral parameter. This can be done only if one of the two weights lu^ or ujq is zero. In the latter 
case, figure [6] shows how to go from the square to the honeycomb lattice. 

There are two possible choices leading to (jjg = 0: u = or u = <& (and other ones which are 
equivalent up to some reparameterization) . The solution u = is trivial as i?(0) oc Id. It is not 
difficult to see that the solution u = $ is related to the 0[n) model on the honeycomb lattice with 
two physical parameters, namely the weight of a closed loop n = —2 cos 4$ and the fugacity of each 
monomer x — 1/(2 cos <J>). Hence, these two parameters are related by 



1 



^J2±^/2- 



(15) 



n 



'We use this parameterization only in this part, because it is more convenient and somewhat usual for the dilute 
loop model on the square lattice. In the rest of the paper we use n = 2 cos 7, which is more convenient when we deal 
with various critical exponents and conformal field theory. 
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Figure 6: When the weight uje = 0, the dilute loop model on the square lattice is equivalent to the 
0{n) model on the honeycomb lattice. 



where the ± sign appears when we consider the different values of $ which give the same n. Now 
we are ready to state Nienhuis' conjecture [19]. The two critical points of the 0{n) model on the 
honeycomb lattice, corresponding to the dilute and dense phases, are the two integrable points. In 
particular, the critical fugacity Xc at the dilute point is 

(16) 



V2 + V2-: 



as claimed in the first part of this article. Actually this was not the argument first used by Nienhuis 
in ^19j, but he formulated it in ^29j . 



3.2 Blob operators 

We will work with three boundary plaquettes 






The piece of loop on the third plaquette carries a blob. It does not matter whether a loop is marked 
with one or several blobs : once a loop is marked, it is blobbed once and for all, and additional 
blobs do not affect this status. In that sense, the blob operator acting on a loop is a projector. The 
physical effect of the blob on a loop is that it changes its fugacity. A blobbed closed loop is given 
a weight ni instead of n, wher^ 

= (17) 
sm 4k$ 



This parameterization is used only in this part dealing with integrability. In the rest of the paper we use 



sin(ri + l)7 
sin ri 7 



and n = 2 cos 7 with 7 g [0, tt), ri e (0, 77/7). 
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Figure 7: Blobs in the dilute loop model on the square lattice. A blobbed loop is a loop which 
carries at least one blob. A blobbed loop has a fugacity ni, whereas a loop without blob gets a 
fugacity n. 



Note that a blobbed loop can touch the boundary without receiving an additional blob each 
time it touches it. It just needs to get at least one blob. If we sum over the possible configurations 
of boundary plaquettes, the total weight of a blobbed loop which touches k times the boundary is 

ni {(32 + Pz)^ — . Pi. loop which is never blobbed would get a weight n/jf • One can redistribute 
those two weights and make things more symmetric if we consider a linear combination of the two 
kinds of loops (blobbed and not blobbed) which carries a weight {n — ni )(32- We already know what 
such an object is : it is the unblobbed loop, which we marked with a square in part El There is a 
slightly more algebraic way to reformulate this. Since the blob acts as a projector on the loops, it 
is natural to introduce the orthogonal projector (or unblob operator) 






and if the three previous boundary plaquettes come with the weights f3i , /32 and /S^ , it is a simple 
change of basis to switch to the blobbed/unblobbed plaquettes. The unblob operator then gets a 
weight P2 , but the weight of the blob operator is (32 + (33 instead of (32 ■ We will use this in section 
13. 3| when we derive the relation (fT2|) for the integrable weights and on the honeycomb lattice. 



There is a nice mathematical object hidden behind this blob operator. The rich underlying 
algebraic structure of loop models is well-known to be the celebrated Temperley-Lieb algebra [30J . 
The extension to loop models with blobs at the boundary is the so-called blob algebra or one- 
boundary Temperley-Lieb algebra. We do not want to elaborate too much about this object here, 
and we refer the reader to references [HU [32l [71 [33] for detailed discussions of this subject. 
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3.3 Solutions to Sklyanin's reflection equation from blob operators 

Recall that Sklyanin's equation for the /C-matrix can be written 




where the dashed lines stand for the identity (or equivalently a contraction of the two edges joined 
by the dashed Hne). Two if -matrices solutions to Sklyanin's equation have been found by Batchelor 
& Yung [17] for the above dilute loop model without the blob operator {Ps = 0) . These two solutions 
are 

Pi = sin(|$ + u) 

P2 = ±sin(|$-u) (18) 
Ps = 

and 

Pi = cos(|$ + u) 

P2 = ±cos(|$-u) (19) 

Ps = 

The -ftT-matrices found by Batchelor & Yung also satisfy the so-called boundary crossing relation, 

which plays some role in what follows. This relation is also satisfied by the two ff-matrices (|24p 
and ([25|l to appear below. 




We still have to explain how we can interpret these plaquettes in the honeycomb limit. Using 
i?-matrices and if-matrices, one can build the family of Sklyanin's transfer matrices for arbitrary 
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wi , a;2, . . . , Wat (we fix iV = 6 for simplicity) 
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(21) 

which commute with each other for any value of the spectral parameter u. Now we choose lji — 
tJa = = w and 0^2 = W4 = = — w, then half of the i?-matrix appearing in Sklyanin's transfer 
matrix ([2T|) are proportional to the identity. 
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The left ii'-matrix can be contracted with the leftmost remaining i?-matrix using (|20p . We end up 
with the so-called diagonal-to-diagonal transfer matrix (we change u into u/2 for convenience) 



Tin) = 



and now it is easy to go to the honeycomb lattice. The important point for us now is that the 
boundary plaquettes take weights /3i(u/2) and /?2(m/2). When we take the honeycomb limit m = <&, 
the second plaquette can be viewed as two boundary half-monomers on the honeycomb lattice. The 
remaining two half-monomers also live on the boundary of the honeycomb lattice, but they come 
from the square lattice model with the bulk plaquettes. To take this into account, it is not difficult 
to check that the weight y of a boundary monomer has to be related to uji and uj2 by y'^/xc = 
/32($/2)//3i(<J>/2). For the two solutions of Batchelor & Yung, this leads to y'^/xc = ±l/(2cos$) 
or y'^/xc = cos$/ cos 2$, which simplifies to 




y 



(22) 



or 



y^ys = {2- n)-^l\ (23) 

The first point corresponds to the ordinary phase (see part [21) • The second one is conjectured to be 
the exact localization of the special transition on the honeycomb lattice. 



18 



The addition of the blob operator in this dilute loop model allows us to build more solutions to 
Sklyanin's reflection equation. Plugging the three boundary plaquettes in the equations, we flnd 
two solutions with 



Pi = sin((2K+ i)$ - It) sin((2K- i)$ 



u 



P2 = sin((2K+i)$ — uj sin ((2k — i)$ — uj (24) 
f]^ = sin(2M) sin(4K$) 



and 



Pi = cos ((2k + i)$ - It) cos ((2k - i)$ 



u 



P2 = cos((2K+i)$-u)cos((2K-i)$-M) (25) 
= - sin(2u) sin(4K$) 

Taking the honeycomb limit u — ^, the weight of a blobbed monomer on the boundary is related 
to the plaquette weights by yfijx = (/32(<i>/2) + /33($/2)) //3i(<i>/2). For unblobbed monomers one 
has Un-m/^ — /32(^'/2)//3i(<&/2). Using the notations of the flrst part (see equation (HO])), one has 

{_ 2cos$sin(2K+l)$ 
~ sin2rt* (ORa^ 

_ 2cos^sin(2K-l)$ {^Dd.) 
^" ~ sin 2k$ 



m 



2 cos <E> cos(2k+1)<E> 

• cos2k<I> 

2cos^cos(2k-1)$ l.ZOU; 



^" cos 2k* 



We conjecture that these two points are the points AS^ and AS,, so they catch the whole universal 
behaviour of the anisotropic special transition. This flnally leads to the relations (fT2l) . 



4 Boundary conformal field theory of the anisotropic special 
transition 

From now on, we turn to the conformal field theory description of the model with the foregoing 
boundary conditions. The conformal field theory of the dilute 0{n) model in the bulk has been 
known for a long time [20]. In particular. Coulomb gas methods have played a crucial role in this 
development. The Coulomb gas approach is a non-rigorous (but powerful) way of constructing 
arguments which lead to exact results, in particular the critical exponents of correlation functions. 
It is a convenient way to deal with loop models, because loops can be interpreted as level Hues of a 
height model, which is argued to fiow towards a Gaussian free field with some background charge in 
the scaling limit. The Coulomb gas approach to loop models is thus heuristically straightforward. 
The Coulomb gas with boundaries, however, is still not well understood [34^ \3E\ [6j. In [6] we 
mixed Coulomb gas with algebraic arguments related to boundary extensions of the Temperley- 
Lieb algebra to derive some results about a dense loop model {ie the dense phase of the 0{n) 
model) with blobbed boundary conditions. The purpose of this section is to extend this approach 
to the dilute loop model. 

4.1 Coulomb gas framework 

4.1.1 The 0{n) model on an infinite cylinder 

Let us start by the usual Coulomb gas approach to the 0{n) model on an infinite cylinder. This will 
fix the coupling constant g of the Gaussian free field as a function of the loop weight n. Consider 
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a loop configuration on the cylinder, and choose an orientation for each loop. Now turn this non- 
locally interacting model (loops are non-local objects) into a local one. Following an oriented loop 
all along, give a weight e*'*''*"/^^'^) to each infinitesimal piece of loop when it turns da, where a is 
the winding angle of the curve. Depending on its orientation, a contractible loop gets a weight e^'''^ . 
The trace over the orientations gives back an unoriented loop with weight n = e*'' + e"*''. Oriented 
loops are then viewed as level lines of a height field h. The height varies by A/i = ±7r when one 
crosses a level line. Then it is generally argued that this model renormalizes towards a Gaussian 
free field h with action 

S=^ f {dhfd^x. (27) 



47r _ 

However, this does not take properly into account the loops which wrap around the annulus, because 
the integrated winding angle / da is equal to in that case, not ±27r. This problem can be solved 
by adding two charges e^^''/'^'''^ at the ends of the cyHnder. This modifies the scaling dimension of 
the vertex operator e'"'* to 

A.-|{(a + 7A)'-(7A)^}. (28) 

The coupling constant g is then usually fixed by the following argument. We started from a discrete 
model in which the height h was a integer multiple of A/i = tt, so the operator cos2/i should be 
marginal. This requires A2 = 2 or A_2 = 2, so g = 1 ± ^. So far, we have not specified if we were 
working with dense or dilute loops. This is what this undetermined sign takes into account. The 
solution 5 < 1 corresponds to the dense phase of the loop model, and g > 1 to the dilute phase. 
Thus, in what follows we will always use the relation between the coupHng constant g and the loop 
weight n 

g = \ + - n = 2cos7 7e[0,7r). (29) 

TT 

4.1.2 Boundaries in the height model 

The arguments here are copied from [6j . To set up a well-behaved Coulomb gas framework, we put 
the model on an infinite strip. The strip has two boundaries, and the key point is to put different 
boundary conditions on the left and right sides. Each time a loop touches the left boundary it can 
get a black blob (or not). A black blobbed loop has a fugacity 

sinfri -|- 1)7 / TT , , , 

ni = ^ ri e (0, -) 30 

smri7 7 

instead of n — 20037. On the right side, loops can get a white blob (or not). A white loop has a 
fugacity 

sin(r2 -f 1)7 , TT, , , 

n2-^7^ r2e(0,-). 31 

smr27 7 

Such a configuration is shown in figure [8l The first step is now to map this model on a height 
model, by choosing orientations for the loops. The question is : what do we do with the blobs? 
Remember that blobs are projectors so if they act on the orientation of the loops, they must be 
projectors on some linear combination of the two possible orientations for one loop. In particular, 
a blob cannot act on the orientation by simply adding a phase shift. Writing that the blob has to 
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Figure 8: Mapping of the loop model with the blobs onto an oriented loop model, and a corre- 
sponding height configuration h/i:. 



be a projector, and using ([30|) — (|3T|) we find that it can be represented as follows 



(32a) 
(32b) 



It is an easy exercise to check that, tracing over the two possible orientations, a blobbed loop takes 
a total weight m (or n2). This works with the convention that a loop has a weight e~*^ when it 
is clockwise oriented. The parameters (pi and ip2 are still free. Actually our problem has a global 
gauge invariance which can be fixed by setting, for example + (/52 = 0. The difference tpi — (f2 is, 
however, very important. It appears when one computes the weight ni2 of a loop which carries the 
two blobs (black and white). Four terms contribute to ni2, because each half-loop can have two 
orientations. This is possible only because the two blobs ((32|) are used here. One gets 



o 



4sinri7sin r27 



-i(ri+r2 + l)7 



„i(ri+r2 + l)7 



-i(i^i-¥'2) 



_ pi(¥'l-¥'2) 



o 



which gives 



sm 



' r-i+r2 + l+ri2 



ni2 



7) sm 



' r-i+r2+l — ri2 . 



smri7sm r27 



(33) 



(34) 
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where we defined ri27 — •^i — for later convenience. 
4.1.3 Flow towards the Gaussian free field 

We have translated the blobs into boundary vertices acting on the orientation of the loops (figure 
[HI). The height model obtained in this way involves complicated local weights on the boundaries 
given by the relations ([32)1 . Fortunately only a few of these weights will be relevant when we go 
to the continuum limit. Our argument for that goes as follows. At the position of a blob, the 
orientation of a loop can be conserved or not (see figure [8]). When it is conserved, relations (|32p 
tell us such configurations come with weights proportional to e^'''^''' or e'^'^^^'^ . How do we take 
these into account in the continuum limit? Our guess is that we do not need to take care of them 
because they do not contribute to the continuum limit. Following [i6j, we argue that since this 
corresponds to the diagonal action of the blobs on the orientations, it can be viewed equivalently as 
a field living on the boundary. Such a boundary condition (with a coupling to a field living on the 
boundary) is expected to renormalize towards a fixed boundary condition under RG, independently 
of the microscopic details about the interaction or the strength of the field. A numerical check 
of the fact that ri and r2 do not appear in the continuum limit will be discussed below (figure 
[T0|) . On the contrary, the parameter ri2 = [^pi — ^P2)/l is important. This parameter appears 
with pairs of vertices (one on each boundary) which do not conserve the orientation of the loops 
(figure [8]). A weight e^^^^'^'^ must be given to each pair of half- loops going from one boundary to 
another. Then we expect that we end up with a Gaussian free field on the strip with action (f27|) . 
with Neumann boundary conditions on both sides, where we have to take properly into account 
these factors e**'"^^'*'. This is the crucial point, which allows us to identify the operator changing 
from one blobbed boundary condition to another. 

Let us be slightly more quantitative by looking at our loop model with blobs on an annulus 
instead of the strip. We want to compute the partition function of this model using the above 
formalism. On the annulus, there can be p pairs of half-loops going from the left boundary to the 
right boundary. This introduces a defect A/i = 27rp when one turns around the annulus. In other 
words, the field h{x,y) can be written as 

h{x,y)^hix,y) + 2npy/T (35) 

where h{x,y + T) — h{x,y) and dxh{x — 0,y) — dxh{x — L,y) — 0. The integration over 
the fiuctuations of h gives the usual factor Zq = q"'^/'^^ / P{q), with q — e^^^^^ and P(q) — 
X\k>\ (f ~ 9'^) • Note that this is just the partition function of the Gaussian free field with Neumann 
boundary conditions on both sides. The other contribution comes from the height defects 27rpj//T, 
which must be counted with a phase e*^''"'*^. The partition function i^o is theiu 

KqO^Z^Y^ g*pri27g-(s/4^)p'(2VT)"(LT) = ^ ^^W^^l 9 1 r)v^ (36) 

where t = T/L. We use the Poisson formula — > J2n I ^P^~^^^"^ ^^"^ perform the integration 
over p. This gives 

Ko oc Zo(T/g)i/2 J2 e-^^-I^B){r,,',i.-2-.)\ (37) 

riGZ 

'^This is actually not really a partition function for the loop model, but rather a Virasoro character. Hence we 
call it Kq and not Z. We will come back to this in section 14.4.21 
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From general arguments we know that such a partition function on a very long annulus (T ^ L) 
should behave as Kq ~ ^o-cju ^j^gj-g y^^ jg ^j^g lowest exponent in the spectrum of the Hamiltonian 
Lq. We use this to normalize properly our result ((37)) . Using the expression of the central charge 



c 1 - 6 
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and Kac' parameterization for the exponents 

{gr - sf -{g- if 



the relation (l37l) becomes 



hr s — 



4ff 



c/24 



nez 



(38) 



(39) 



(40) 



which should be interpreted in boundary conformal field theory (BCFT, see [Ml E?]) as follows. 
First, we interpret our Gaussian free field as a 1 + 1-quantum field theory quantized on a segment 
of length L. In this language, the Hamiltonian H is the generator of translations in the T-direction 
on the annulus or on the strip. The strip can be mapped by the conformal transformation z i— > 
exp(— -^2:) onto the half-plane (figure [9]), where H is related to the usual Virasoro generator Lq 
(the dilatation operator) by the well-known formula H — (Lq — ^) . The spectrum of H and the 
spectrum of Lq are then related. But Lq is a Virasoro generator: it acts on some (highest-weight) 
representations of the Virasoro algebra. In conformal field theory, the highest-weight representations 
are the primary operators. Where are these operators in the half-plane? They must be situated at 
the origin, in order for the boundary conditions to change between the positive and negative real 
axes. Such operators are boundary- condition- changing- operators (B.C.C) and were introduced by 
Cardy in [37]. Computing the partition function ([40]) on an annulus of size L x T is equivalent to 




exp [ — ^z) 




Figure 9: The infinite strip is conformally equivalent to the half-plane. When the boundary condi- 
tion is different on both sides of the strip, a B.C.C. operator is sitting at the point 0. 



taking the trace of the evolution operator e . Recalling that we defined q — e ^^/^, such a trace 
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should have the form tr {e "^^} = (l/^(<z)) 'J2h q^"^''^^^ where {ha} are all the primary operators 
appearing in the spectrum of Lq. The factor 1/P{q) comes from the trace over the descendants of 
each primary operator. This is indeed the form of formula ((40)) . 
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Figure 10: The leading exponent ho in the sector without strings is written ho = (C^ — l) ^^^^-^ . 
We plot C vs. ri (a) and ru (b). Numerically, it is clear here that $ does not depend on ri at all, 
and that ^ = ri2, as predicted from the Coulomb gas arguments. 



We have thus identified the spectrum of primary operators which are necessary to change the 
boundary condition from the black blobs on the left to the white blobs on the right. These primary 
operators depend on the loop weights ni, n2 and ni2 only through the parameter ri2 and the 
relations ((30l) . (|3ip . (l34)) . according to the above argument. We have checked these results numer- 
ically by transfer-matrix diagonalization of the dilute loop model on the honeycomb lattice. The 
free energy per unit site /l is computed for successive widths L — 8, 9, . . . , 12. We then use the 
well-known relation 

/. = /bu. + + ^Jh^^g^ + o (i, ) (41) 

up to order O (-^) to extract the lowest exponent ho appearing in the spectrum of Lq. The result 
predicted above ([40]) is that ho = ^•ri2,ri2 for ri2 € [0,77/7). Results are shown in figure fTOl for 
n — \/2 (so ^ — J and c = ^) and for anisotropic special boundary conditions AS, and ASn on 
the left and right sides. They are in very good agreement with the Coulomb gas prediction. 

4.2 Boundary operators : B.C.C versus string-creating operators 

4.2.1 From ordinary to anisotropic special boundary conditions 

In section 133] we computed the spectrum of the B.C.C operator changing from one type of blobs 
to another (formula (|40|) ). and we saw that the weight nu given to loops touching both boundaries 
was a crucial ingredient. We would Hke to use this information to understand a simpler B.C.C 
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operator : the operator going from the ordinary boundary condition to the anisotropic special one. 
To do this, we must give to the (white) blobbed loop the weight 712 = n, and to the loops carrying 
the two blobs the weight ni2 = ni. In that case there are no closed white unblobbed loops, because 
they would get a weight n2 — n — 0. If a loop is blobbed it does not get an additional weight each 
time it gets a new blob on the boundary because Wo ^ I (see relation (fT2|) for the point ASa) ■ Thus 
we identify the ordinary b.c with the point ASa when n2 = n. Now recall the parameterizations 
((3T1) . (|34p to see that r2 = 1 and ri2 = ri. Then the leading exponent appearing in (|40)) is hr^^n- 
Our guess is that the primary operator $ri,ri is the B.C.C operator from ordinary to anisotropic 
special. Note that, since in general ri does not need to be an integer, <i>ri,ri is the highest-weight 
state of a generic Verma module, ie without null states. 

However there is something we did not specify yet: what anisotropic special b.c we are talking 
about. Is it AS^ or AS^l To answer this, note that all the different quantities in our model are 
supposed to vary continuously when we vary ni. Taking the limit ni 71 we see that the B.C.C 
operator is the identity so the boundary condition is actually the same on both sides of the B.C.C 
operator. AS^ is not the ordinary b.c. in that case (actually it is the isotropic special one), whereas 
AS, fits. Hence our conclusion is that for arbitrary ni the B.C.C operator {AS^/Ord) is $ri,ri- 

The operator changing from ordinary to the other anisotropic special b.c can be identified by 
symmetry. Indeed, note that the model ^ is invariant under the change ni ^ n — ni and • ^ ■ in 
all the formulas. For example, the points AS» and AS^ are exchanged. Thus the B.C.C operator we 
are looking for is the same as the previous one if we replace ni by n — m . This is the same as taking 
7/7r — ri instead of ri. Using Kac' parameterization ([39|) we find that ^7/7r-ri,7/ir-ri = ^ri,ri+i so 
we conclude that the B.C.C operator {AS^/Ord) is <i>ri,ri+i- 

4.2.2 Strings 

In a loop model, a very natural boundary operator is the operator creating a piece of loop (a string) 
at a point on the boundary. A string cannot stop in the the bulk or at a boundary point, therefore it 
must propagate to infinity. The string-creating operator is known to be 4*2, 1 with scaHng dimension 
h2,i [SI EH El. It has a null stat(ll at level 2 

(^J:-2-^i-i) |*2,i> =0 (42) 

the link between this operator, the relation (|42|) and the Schramm-Loewner evolution SLE^, n — - 
is presented in [39|. One can generate more than one string at a boundary point. However, when two 
strings are generated, they can be contracted into a closed loop, without propagating to infinity. 
On the contrary, one can ask that they both reach infinity, but this corresponds to a different 
representation of the Virasoro algebra. This can be formulated in terms of fusion as 

^'2,1 ® f ^'24 = © *3,1 (43) 

*We use the convention that a primary operator ^r,s has no null state (and hence corresponds to a generic Verma 
module) whereas '^r,s with r, s e N has a null state at level rs (hence a submodule of the generic Verma module is 
quotiented out). 
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where the operator ^i^i is just the usual vacuum of the theory (this is the case when the two strings 
are contracted) and '^3^1 creates two strings, conditioned not to annihilate with each other. Note 
that ^3.1 has a null state at level 3. Let us give a more pictorial view of this: 




*2,1 *2,1 ^'3,1 



The operator creating L strings can be found by induction, iterating the fusion with ^'2,1- It is of 
course (with a null state at level 1 + L). 

What become the strings when we put the anisotropic special b.c AS, or AS, on the left part of 
the boundary? The result has to be interpreted in terms of the fusion of the string-creating operator 
5*2,1 and the B.C.C operator ^ri,ri+i (or $ri,ri in the case of ASm)- Since 5*2, 1 has a null-state 
at level two, the fusion can be computed by writing down the differential equation satisfied by the 
three-point function {^ri,ri+iizi)^siz2)'^2.i{z3)) to identify the possible scaling dimensions of ^s- 
The result is well-known 

*ri,n + l «> / ^'2,1 = *ri + l,ri + l © $ri-l,ri + l (44) 

The two terms on the right correspond to the two possible states of the string generated by "^2,1, 
which can now be either blobbed or unblobbed. In the case of the b.c AS,, the weight given to a 
blobbed monomer is bigger than the one of an unblobbed monomer. Then the scaling dimension of 
the operator creating a blobbed string should be smaller than the one of the operator creating an 
unblobbed string in that case. In general, hr-^+i^n+i < ^ri-i,ri+i so we get the following pictorial 
view 




One can add more strings by iterating the fusion procedure. However, note that only the leftmost 
string can be blobbed or unblobbed, because the other ones never touch the boundary. The operator 
creating L strings with AS, b.c on the left and with the leftmost string blobbed is then ^n+L.n+i- 
When the leftmost string is unblobbed, it is <i>ri-L,ri-i-i- In the case of AS, b.c on the left, the 
operator is ^n+L.ri when the leftmost string is blobbed, and ^ri-L,ri when it is unblobbed. Note 
that again, the rii ^ n — ni duality shows up, because all these operators are exchanged when one 
changes all the • into ■ and n into ir/j — ri. 

4.2.3 Virasoro characters 

So far, we have identified the B.C.C operators going from ordinary to one of the anisotropic special 
boundary conditions, as well as their fusion with the string-creator operator '^2.1- Taking the 
trace of qLo-c/24 ^^er the Verma modules of these operators, we obtain the corresponding Virasoro 
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characters. L is the number of strings and we indicate the blob status of the leftmost string when 
necessary. K^f is the character without string {L — Q). For completeness we indicate the characters 
for anisotropic special b.c on the left and right side. These are conjectured from a combination of 
results of [6] and from algebraic arguments coming from the two-boundary Temperley-Lieb algebra. 
We do not want to develop this here, and relegate it to [33] . 



Ord/Ord : 



AS^/Ord 



KUq) 



AS,/Ord : 



AS,/ASo 



Ko{q) 
Kl°{q) 
Kl°{q) 
Kl°iq) 
Kl°{q) 








(45) 



(46) 
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1-2T1 



P(q 



l-2n 



P(q) ^r. 
c/24 



~ P{q) l^n>oq 



-ra-l + i.ri 
/l_ri+r2-l + -L,-ri + r2-l-2n 
2-1 + 1,, -i-i-T-2-3-2n 



(48) 



The cases of ASu/ASn, ASm/ASo, ASm/ASn can be deduced from the AS^/ASo case, by the 
duality transformations ni n — ni and/or 712 — "2, as explained previously. 



4.3 Some remarks about fractal dimensions 

The link between scaling dimensions of operators and fractal dimensions of some geometric objects 
is well-known. For example, in the dilute loop model with which we are dealing here, the fractal 
dimension d/ of the loops in the bulk can be computed very easily. It is the RG eigenvalue of the 
2-arm exponent hi,o, which gives the celebrated formula df = 2 — 2hi^o ~ ^ ~^ J^- '^^^ boundary 
operators we have identified above might be used to derive such formulae. The boundary condition 
cannot change the local geometry (and especially the fractal dimension) of an object at a point in 
the bulk, so here the fractal dimension of the loops is not the quantity we are interested in. We 
should rather look at the set of points on the boundary which are visited by a chosen loop. Let us 
call this the contact set of a given loop. As above, the fractal dimension d^^''^'^ of such a set must be 
given by a (boundary) 2-arm exponent, which depends on the boundary condition (the same on the 
left and right sides of the point where the operator is inserted). But we have just identified these 
operators for all the boundary conditions we are interested in. For example, for ordinary boundary 
conditions, one has 



,{Ord) _ 



1 



i3,l 



< 



(49) 
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Figure 11: The contact set (blue) of a loop is the intersection of its curve 7 and the boundary of 
the upper half-plane (the real axis) . Its fractal dimension is related to the scaling dimension of the 
2-arm exponent. 

for < n < 2. This means that, in the scaling limit, the contact set of a given loop is almost surely 
empty: loops never touch the boundary. On the contrary, at the (isotropic) special point, this set 
has a non-trivial scaHng-Hmit with fractal dimension 

= 1-/^3 3 = 1-2^^^ (50) 

When the b.c is AS, there are two kinds of loops: blobbed and unblobbed ones. The contact set of 
both types are different. For a blobbed loop the fractal dimension is 

A"^^'^ (. loop) = 1 - /l2ri + l,2ri + l = 1 - n (n + 1) ~ (51) 

' 9 
while for an unblobbed loop 

df^*^ (. loop) = 1 - /i-2n+i,-2n-3 < (52) 

for < ni < n < 2 so in the scaling limit, an unblobbed loop never touches the boundary. We 
conjecture that these fractal dimensions might be of some relevance for SLE^^ p when 8/3 < k < 4. 
So far (to our knowledge), rigorous results have been obtained only for the contact set of SLE^. 
when 4 < K < 8 gO]. 

4.4 Annulus partition functions and boundary entropies 

We use the above Virasoro characters to derive expressions for the annulus partition functions of 
the 0{n) model with ordinary and anisotropic special boundary conditions. The geometry of the 
annulus plays a role in BCFT similar to the one of the torus in CFT. A partition function can be seen 
either as a periodic evolution of a one-dimensional quantum theory with two boundaries (also called 
the open-string channel) or as the evolution of a periodic theory between two boundary states (this 
is the closed string channel) [37]. So far we have dealt only with the open-string channel, working 
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with B.C.C operators in the half-plane, then deducing the corresponding Virasoro characters. In 
this section we conclude this work by gathering all the pieces of the puzzle to obtain the annulus 
partition functions. Then we go to the closed-string channel and compute the boundary entropies 
of the boundary states 10?^^), I^S",), etc. 

4.4.1 A loop partition function is a Markov trace (not a trace) 

A loop configuration on the annulus can be drawn as a planar rectangular diagram with half-loops 
attached on the top and the bottom of the diagram. When a closed loops appear in the diagram, 
it can be removed if we give a weight n to the diagram. With this rule we get a set of diagrams 
connecting N sites on the bottom to N sites at the top of the diagram. Two diagrams are equivalent 
iff they represent the same connectivities between the sites. For example 




If we define the product of two diagrams as a concatenation, putting one diagram above the 
other and counting the closed loops with a weight n, then this defines the Temperley-Lieb algebra 
TL N {n) . The generic representation theory of this algebra is very simpl^. A generic irreducible 
representation is given by the set of half-diagrams connecting N points. The diagrams of the TL 
algebra act on these by concatenation. For example the set 

vo = { W W , ^} 

is a representation of the algebra TLi{n). Other irreducible representations are obtained when one 
introduces strings. A pair of strings cannot be contracted: a TL diagram which makes a connection 
between two strings takes a weight zero. With this rule, one obtains two additional representations 
for TL4(n) 

V, = {\J I I , I U I , I I U } and V4 = { I III} 

Now we are ready to the introduce the Markov trace. The Markov trace of a diagram is the number 
obtained when the top and the bottom of the diagram are identified and each closed loop is given 
a weight n. For example, the Markov trace of the above diagram M is Tr M = n^. The Markov 
trace extends to the combinations of TL diagrams by linearity. Note that it is also invariant under 
cyclic permutation of the diagrams: for two diagrams A and B which are multiplied with the 
concatenation rule, Tr AB — Tr BA. Note that the Markov trace is exactly the object we need to 
compute a partition function of a loop model on an annulus. The configurations of the loop model 
are given by the TL diagrams with suitable weights, and the geometry of the annulus is obtained 
when one identify the top and the bottom of each diagram, with a weight given by the Markov 
trace. 

^Generic means when n is not of the form 2 cos (^) with m an integer. In the non-generic case most of the string 
representations become indecomposable. The irreducible representations obtained by taking quotients of the latter 
ones are related to the RSOS models 41_. This comes out in the Markov trace because some coefBcients are zero in 
that case. 
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The Markov trace of a diagram is very different from the usual trace over a representation of the 
TL algebra. However, there is a well-known relation between these objects 

TrA^^PLtrv^A (53) 

L>0 

where the Vl are Chebyshev polynomials of the second kind I?o = 1, T^i = I?2 = ri^ — 1, 
2?3 = n{-n? — 2), etc. In full generality = '^^^^ ^ ~ 2cos7. Of course, we use the 

convention that a representation Vl of TLiq{n) is empty if L and N do not have the same parity, 
or if L > N . Again consider the example of the diagram drawn above: it is easy to check that 
tr^gM — n^, try^ = and tr^^ — and all the other traces are zero by definition, so using the 
formula ([53)1 one finds Tr M = n^ as expected. 



It is clear now why one has to be interested in the relation ([53)1 in view of the preceding sections. 
Indeed, we have identified the natural candidates for the different string representations of the 
TL{n) algebra in the scaling limit. These are the Verma modules of the string-creating operators 
'^i+L,i- Thus in the scaling limit, the trace over the different TL(n) representations of the evolution 
operator (on the lattice this is some power of the transfer matrix), when properly renormalized, 
converges to the trace over a Verma module of the BCFT evolution operator q^a-c/2A^ These traces 
are the Virasoro characters (|45)) . Thus we get the partition function of the 0{n) model on an 
annulus with ordinary b.c on both sides for free 

Zord/Ord = Y.'^LKL{q). (54) 

L>0 

Interestingly enough, the fact that the coefficients Vl are Chebyshev polynomials of the second 
kind has a nice interpretation in terms of fusion. Indeed, we have 

VlVi^Vl+i+Vl-i (55) 

and this looks like the fusion rule (see (|43|) for the case L = 1) of the operators and 4*2, i. 

Another way to say this is that the coefficients Vl, also named quantum dimensions, satisfy the 
fusion algebra of '^i+l.i- This gives a simple and powerful way of computing the coefficients Vl- 



So far, we have not talked about the blobs in the present section. Actually the generaHzation of 
formula ([53]) to the case with blobs is straightforward. The whole framework of the Temperley-Lieb 
algebra can be generalized to a diagrammatic algebra where some of the loops can carry blobs 



niTL 




This algebraic structure has first appeared in ^T] under the name of hloh algebra, and it was 
sometimes renamed into one-boundary Temperley-Lieb algebra or IBTL since [32l [7l [6j l4j [5] . This 
algebra has now two parameters n and ni. The generic representation theory follows the one of 
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Figure 12: Open-string channel (on the left) and closed-string channel (on the right). 



the TL algebra. The representations are the sets of half-diagrams connecting N points without 
crossing. Some of the half-loops can now carry blobs. In the case of representations with strings, 
there are two inequivalent representations, depending on whether the leftmost string is blobbed 
or unblobbed. There exists a Markov trace over this algebra, and an equivalent of formula (|53l) 
involving all the string representations (blobbed and unblobbed). The corresponding coefficients 
and I?" turn out to satisfy the fusion rule (|44|) . or in full generality |4l[5] 

VlV, = + Pt-i = n+i + (56) 

and Vq = l,!^* = ni,T>^ = n— ni. This can also be written as = '^'"in^^^'*'^ and I?" = ^^^^^777^^- 
Again, with this result at hand, we can write down partition functions very easily. For example, the 
partition function of the 0{n) model on an annulus with ASu/Ord boundary conditions is simply 
given by the characters (f46|) and 

ZAs,/Ord - Ko{q) + J2 T^lKliq) + ^ P^/f^(g) (57) 

L>1 L>1 

In the case of two boundaries, there is also an equivalent of formula ([53)1 . this time involving all 
the string representation with the different blob status of the leftmost and rightmost strings. More 
details about this, as well as the expressions for the quantum dimensions 2?"°, etc. can be found 
in 0[5]. 



4.4.2 Results for the boundary entropies 

We have just seen that it is now very easy to write down the partition functions of the 0{n) model 
for the various boundary conditions we are considering. We use these partition functions to compute 
the boundary entropies [42] of the corresponding boundary states. Recall that in BCFT, a partition 
function on the annulus can be viewed either in the open-string channel or in the closed channel [37] 
(figure [T2I) . So far we have worked only in the open-string channel, where our partition function are 
Markov traces of the evolution operator q^o-c/24 ^-^ _ e"^"^). We go to the closed-string channel by 

performing Poisson resummations / dLe^^'^^^ on the partition functions. The result is of 

L p 

the form Z — {a\ qLo+Lo-c/12 where q = e^'^'^l'^ is the modular parameter conjugated to g, and 
I a), 1/3) are the boundary states which implement the boundary conditions a and on both sides 
of the annulus. In the limit of a very long cyHnder (g ^ 0), only the lowest exponent contributes in 
the closed-string channel. Here the lowest exponent is zero, so Z ~ («|0) where |0) is 

the ground state of Lq = ^0 — c/12. The partition function factorizes into contribution from each 
boundary, and a bulk part which is trivially 
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The g-factors ga — (a|0) and gp — (/3|0) count the degeneracy of the ground state in the 
boundary states \a) and [42j. The free energy / — — \ogZ gets a contribution /i.b = — log^o, 
from the left boundary, and fr.b = — log^a from the right one. These quantities are called the 
boundary entropies of the boundary conditions a and (3. The g-factors and the boundary entropies 
are of interest in the study of BCFT. We will use them in section [5] to organize the boundary RG 
flows between the different boundary conditions of the 0{n) model. 



Ordinary boundary condition on both sides: the partition function is l|54p 



Zord/ord-^^}^^ sin7 



and its modular transform is 



/9-;;-i/i2 sin(^-^27rn) 

^Ord/Ord - \ 2^ 



(2p+(,-l)" 



9 Pm ^ sin 7 

which gives the boundary entropy Sord — ~ log gord where 



(Isotropic) Special/Ordinary b.c: The operator going from ordinary to (isotropic) special b.c 
is known to be ^'1,2, so the Virasoro character with L strings is 

-c/24 

ifl + L,2 = ^(9''^ + -^-'z"^+-'-) 

and the partition function is then 

^Sp Ord — ^ / , ■■ [Q - q ) 

P{q) ^0 sm7 ' 



the modular transform is 



/9- ~-i/i2 sin (2^ + ^2Trp] , 

Zsp/Ord^J--^^ — 



q P(q^) ^ sin 7 

which gives the boundary entropy Ssp — — log gsp with 



, ^y/^ sin(27/.g) 
^y/ (sm7sm(7/g)) ' 
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Anisotropic special / Ordinary b.c with the anisotropic special b.c AS^ on one side of the 
annulus and ordinary b.c on the other side, we obtain from (|46P and ((57)) 



q ^/^^ ^ sin(ri +L)7 , 

and its modular transform 

9~-i/i2 sin ( in + 1)7/0 + ^-'•itg'^) 27rp') 



- 

pez 



\/ g p(g^) smri7 



so we have the boundary entropy Sas^ = — log gAS^ , with 

^2V/ Sin((ri + l)7/g) f sinj 

9AS, = - ■■ ■ / / \ ■ 60 

\gj smri7 Vsm(7/.g)y 

The boundary entropy Sasm = ^ ^'^S 5^5, can be deduced from this one by the duality transforma- 
tion ri — > 7r/7 — ri. We get 

9AS.^(-X'''^^^^^(-^X'\ (61) 



g J sinri7 \sin(7/(7) / 



5 Boundary RG flows 

The BCFT machinery developed in section [4] is now used to explain some features of the phase 
diagram shown in figure HI We also discuss the link with known results of integrable field theory 
(IFT) and the thermodynamic Bethe Ansatz (TEA). 

5.1 Affleck &: Ludwig's "(^-theorem" and consistency of the phase dia- 
gram 

In f42], it was argued that the g-factor (see section l4.4.2p is always decreasing under renormalization 
from a less stable to a more stable boundary condition. The g-factor plays a similar role in BCFT 
to the one of the central charge in the bulk CFT [24]. Again, this is true only when the theory 
is unitary, which is not the case in general for the 0{n) model. However, it turns out that this 
still makes sense here and that we can organize the boundary RG flows according to the boundary 
entropies of the different boundary conditions we have considered so far for the 0{n) model. The 
g-factors l|58)) . (|59p . ((60)) . ((6T|) are plotted in flgure [TSl for n = 1. We see that the RG flows between 
the different b.c predicted by Affleck & Ludwig's argument are in agreement with the phase diagram 
Hlin the region < ni < n. The most unstable b.c is the (isotropic) special one. The two anisotropic 
special b.c are more stable than the special one, but less than the ordinary b.c. 

Note that outside the region < ni < n, the ^-factors are organized differently, so the phase 
diagram[4]is probably wrong in this case. Therefore, in this article we restrict our discussion to the 
"physical" region < ni < n. 
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Figure 13: g-factors of the ordinary, special and anisotropic special boundary conditions for the 
0{n) model. The results here are shown for n — 1. In the "physical" region < ni < n, the 
^-factors suggest that the most stable b.c is the ordinary one, and the most unstable is the special 
one. This is in agreement with the phase diagram shown in figure |4l 

5.2 Boundary perturbing operators : some guesswork 

We want to understand why the ordinary b.c is more stable than the special one when one perturbs 
the theory in the anisotropic direction. This can be done within the framework of perturbed BCFT. 
The idea is that, at an RG fixed point, the theory is described (formally) by an action Sbcft- One 
can perturb this theory by a boundary operator $6 with scaling dimension /i^ 



the eigenvalue of the coupling A is then y\ = 1 — hb under the RG fiow. If one identifies the operator 
$6 properly, one can then say if the perturbation is relevant {y\ > 0) or not. 

In the case of an anisotropic perturbation around an isotropic b.c, the perturbing coupHng is 
A oc A where A = — is called the anisotropic coupling. Now, what follows is very close to 
what we said about the fractal dimension of the contact set (see previous part) for special and 
ordinary b.c. Indeed, an operator perturbing in the anisotropic direction should create additional 
loops on the boundary. We already know what operator creates a piece of loop at a boundary 
point: it is the 2-arms operator. In the case of ordinary b.c, it is <i>3,i, which is irrelevant when 

< n < 2. But for special b.c, it is $3,3 and this is relevant, with eigenvalue j/a = 1 — g ■ 
We can reformulate this in terms of the fractal dimensions discussed above. For ordinary b.c, loops 




boundary 



(62) 



34 



in the scaling limit almost surely never touch the boundary. Hence, a very small perturbation at 
the boundary could hardly affect their critical behaviour. At the special point, however, the loops 
can touch the boundary in a non-trivial way (on a contact set with a non-trivial scaHng limit) . For 
such a boundary condition, a small perturbation at the boundary can affect the critical behaviour 
in a more crucial way. 

There is another perturbation of interest around the (isotropic) special point: the perturbation 
in the isotropic direction. The coupHng in that direction is the weight of a boundary monomer y 
(recall it is also the coupHng between two boundary 0{n) spins, see relation ^). In that case the 
perturbing operator has been known for a long time [43j to be <&i,3, with eigenvalue Uy — 2^^^^-. 

5.3 Shape of the phase diagram around the isotropic special point 

In this section we explain why the phase diagram shown in figure [4] has a cusp-like shape around 
the special b.c. To do this, we apply finite-size scaling arguments [21} 144) . Consider the boundary 
free energy /fe(y, A) as a function of the isotropic coupling y and the anisotropic coupling A around 
the special b.c y = ys and A = 0. The singular part of ft has the scaling form 

/m(2/,A) = 1vI.(^^,_^^ (63) 

where ^f, is the correlation length along the boundary (the bulk correlation length is infinite because 
X = Xc). Restricting to y < ys, one can apply renormalization group transformations until we reach 
{ys - y)ib^'' = 1 and then 

There is a similar function for y > ys but we already know that in that region one goes from one 
extraordinary phase to another when A = and that this is a first order transition (see section [2]) . 
In the case y < ys we should cross the anisotropic special transition line for some A+(y) > and 
A~(?/) < 0. At these two points fb^s must be singular, whereas it is non-singular in the rest of the 
domain. We have then A* (y) = C^{ys ~ y)^^"^, with = yy/y/^- According to the previous section 
(f) = ![~^^'|^ < 1 for < n < 2, so the anisotropic special transition must have the cusp-Hke shape 
drawn in figure IH 

We have checked the value of the exponent 4> numerically by transfer-matrix diagonalization. We 
compute the ground state of the transfer matrix on a strip of width N . We compute the free energy 
per site /at, which should behave as (|4T|) . Then we extract the effective central charge CcS = c— 24ft,o 
of the theory. Once again we introduce the quantity C related to the ground-state exponent /iq by 
ha = (C'^ — 1)^^^^. We expect a scaHng behaviour of C, of the form (see figure [HI) 

It must be possible to relate the exponents 5y and (5a to the usual exponents of the boundary 
correlation lengths (parallel and orthogonal), and to yy and yA- However, because of the mixture of 
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Figure 14: The leading exponent Hq in the sector without string is written ho = (C^ — l) ^ ■ 

We plot (vs. y and A around the special point. Here we have 71 = \/2, so the central charge is 
7/10 and g — 5/4. At the special point ho — /ii,2 (so here = |5 x 1 — 4 x 2| = 3). Extraordinary 
b.c give ho = so C = 0, and extraordinary b.c ho = /i2,i so C = |5 x 2 — 4 x 1| = 6. The results 
here are in good agreement with 5y = 0.49 ± 0.02 and Sa — 1-1 ± 0.05. 



boundary and finite-size effects [HIH], this would be compHcated. Fortunately we do not have to 
go through this, because if we simply take the limit TV ^ 00, for y < ys we arrive at the following 
scaling form 

where O must now be a piecewise constant function. Then we see that = so even if we are not 
able to relate yy to Sy and yA to Sa independently, we can still compute there ratio (j) numerically. 
We find 
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g 
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Sa 




(f) (theor.) 





3/2 


0.87 ±0.05 


1.00 ±0.05 


0.87±0.1 


1 


1/2 


4/3 


0.63 ±0.03 


1.07 ±0.05 


0.59 ±0.08 


3/5 = 0.6 


7/10 


5/4 


0.49 ±0.02 


1.10 ±0.05 


0.45 ±0.07 


4/9 = 0.444 


8/10 


6/5 


0.41 ±0.02 


1.25 ±0.05 


0.33 ±0.07 


5/14 = 0.357 



The results are in good agreement with the relation (j) — ■ Note however that the limits n ^ 

and n — > 2 give results which are less accurate than for n = 1 or n = \/2. 
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5.4 Integrable flows and TBA 



So far, we have not discussed the perturbation around a point on the anisotropic special Une. Note 
that the whole line is attracted by one of the two points AS, or AS,, so it is sufficient to study 
only the transition at one of these points. In particular, the crossing behaviour along the whole line 
must be described by the perturbation around AS, or AS, in the unstable direction. This is what 
we study in this section, by means of the Thermodynamic Bethe Ansatz (TBA). 

The TBA system for the 0{n) model in the bulk and in the rational case n — 2cos(^) was 
derived in [45]. In the massless case of interest here, it becomes 

where j = 1, . . . , m — 2 and eo — em-i = oo. 



In [46] Fendley derived the effect of the introduction of a boundary S'-matrix of arbitrary spin. 
It turns out to add a term in the free energy of the system 

= £ cosh(.-ilg(r/r.)) ^"K^ + -"^^^^^0 

where, in the language of the Kondo problem, S is the spin of the magnetic impurity and Tk is the 
Kondo temperature. 



We analyze this TBA system (HH-lISHl) following [43]. First, let us focus on the UV limit T > Tk- 
The leading contribution to the boundary free energy ([68|) comes from the region 9 close to infinity. 
Introducing Xj = e~'^^^°°^ and taking the Hmit 6* — > oo in ([GT]) . one gets 

= (l + xj--i)(l + a;j+i) (69) 

for j = 1, . . . , TO — 2, xo = Xm-i = 0. This yields 

l + .,= f""".(^'+;^/("^+^^V. (70) 



sinTr/ (m + 1 

The leading contribution to the boundary entropy is then 



AUV) 

■/boundary ^ /i , ^ , /-/^ / ^ \ ^ 

jr^ = -^log{l + x2s) + Ol — ]. (71) 



Now, let us turn to the IR limit T <C Tk- In that case the leading contribution to (|68|) comes 
from the region 9 — > — oo. Introducing yi — e~'^^^~°°\ the 9 — oo limit of system l|67|) gives 

yi^O y| = (l + y,-i)(l+2/,+i) (72) 

with j = 2, . . . , TO — 2, Um-i — 0. This is solved by 

/ simr j/m\ ^ , , 

\ smn/m J 
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Now the leading contribution to the boundary entropy is 

f^"^'> 1 / T \ 

•^boundary ^ fi , \ I \ /r^A\ 

— 7^ = -7^^od^ + y2s) + o\ — \. (74) 

In this TBA framework, the g-factor is related to the boundary free energy hy g = exp (— /boundary /T) 
times some overall constant which is independent of the scale T/Tk- Thus the TBA system ([67)l -(|68 )) 
leads to 

9^ _ (1 + X2g)^/^ _ sin7r(25' + l)/(m + 1) sin7r/m 

giR ~ {l + y^sY/-^ ~ sin7r/(m + l)sin7r(2S')/m ' ^ ' 

Comparing this to (|58p and ([60)) . we see that with ri = 25* one has 

4^ = ^ (76) 

so the TBA system above is describing the flow from the anisotropic special b.c with ri = 25* to 
the ordinary b.c. The case ri = 1 is actually the case of the flow from the (isotropic) special b.c to 
the ordinary one, and this was studied in |43j. One can also see here that there is no chance that 
we get in the end the flow from (isotropic) special to anisotropic special b.c, so this flow is probably 
not integrable. 

The TBA system exposed here allows us to compute the dimension of the perturbing operator 
around the anisotropic special b.c. Again, this can be done following [43]. Consider the UV limit 
T ^ Tk- Then the TBA system implies the periodicity of the pseudo-energies ej{9 + (m + l)i7r) = 
ej{9), so close to 6* = oo we can expand 



log(l + e-^^(«)) = 4'^ (e-^e/^^+i)) ' (77) 



A;=0 

Plugging this into ((67| . one can show that the term fc = 1 is zero. Then in the UV limit the 
boundary free energy can be expanded as 

f(UV) oo s 2k/(m+l) 

■/boundary ^1/1, \ , \ AUV) /„„x 

so the scaling dimension of the energy operator at the anisotropic special transition is always 
= 1 — Thus the operator perturbing in the unstable direction around the anisotropic 

special transition is $1,3. 



6 Open boundary conditions 

In this section we focus on "open" boundary conditions for the 0(n) model. At flrst glance, these 
conditions are different from the ones we considered so far. However, it turns out that these 
conditions are related to the anisotropic special b.c in the particular case ni = 1. It is then possible 
to use the foregoing formalism to study these b.c in the BCFT framework. As an application, we 
derive a crossing formula for the Ising spins clusters on an annulus in section 16.41 
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6.1 Cardy's open boundary conditions, boundary fields in the spin 0(n) 
model 

In an unpublished note [47], Cardy introduced the b.c for the 0{n) model shown in figureflU which 
we call open b.c. The boundary of the honeycomb lattice is sHghtly different from the one we were 
looking at so far. This is important since it is now possible to have half-loops attached to the 
boundary. A half-loop gets a Boltzmann weight vi. 



a 
o 




Figure 15: Open boundary conditions in the 0{n) model. The closed loops in the bulk still have a 
fugacity n, but the half-loops have a fugacity i^i. 

In more physical terms, one can think about these open boundary conditions as a magnetic field 
Bi at the surface, which couples to the boundary spins. 



^ = Tr 



boundary bulk <ij> 



(79) 



where Bi is a boundary field. Proceeding as usual, we expand this expression and take the trace 
over independent terms to get the familiar loop expansion, with half-loops ending on the left and 
right boundaries. These loops get weights 



1^1 



B, 



(80) 



6.2 Relation with the anisotropic special b.c when ni = 1 

Open boundary conditions correspond to a surface magnetic field coupling to the boundary spins. 
We expect such a boundary condition to fiow towards fixed boundary conditions in the 0{n) model 



as soon as 



Bi 



> 0. This leads to the simple result that, in the scaHng limit, the b.c condition 
should not depend on the precise value of vi as soon as > 0. It is easy to check that fix b.c 
for the 0{n) model are exactly the open on^l^ with vi — 1. Thus the whole critical behaviour for 



^"Although it is clear from this argument that the open b.c is simply the fixed b.c in the scaling limit of the 0{n) 
model, we keep calling it "open" for the consistency of our terminology. 
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1^1 > is completely caught by the theory at = 1. 




Figure 16: In the special case ni = 1, there are no unblobbed loops and the weight of the blobs 
compensates exactly the weight of the boundary monomers. A loop configuration on the lattice (a) 
is then equivalent to a configuration on the lattice (6) with open b.c when ui = 1. Since all the 
open boundary conditions with vi > Q are expected to renormalize to the same b.c, they must all 
be described by the theory with ni = 1. 

Now consider the anisotropic special b.c AS^ in the case ni = 1. Note that the formulae (|12p 
simplify in that case to = and = 0. There are no unblobbed loops because — 0. A 
loop touching several times the boundary gets a weight WuX^ = 1 per pair of boundary monomers 
(see figure [TH). Then we can just erase the first row of spins of the honeycomb lattice, and we end 
up with the lattice which gives rise to the open b.c, here with — 1. Because of the above remark, 
this describes the open b.c for every vi > 0. 

It is easy to apply the results of section |4] to derive new results about the BCFT of the open 
b.c. The constraint ni = 1 fixes ri — 2^^^—!-. Then the B.C.C operator going from ordinary b.c to 
open b.c is ^n.ri+i, which means that its scaling dimension is 

h ^/~,-l , -I = /i^tAt + I tt/t + I (81) 

2 ■ 2 2 ■ 2 

where we used the symmetry of Kac' formula ([39]) . This operator, which should be viewed as the 
B.C.C operator going from Dirichlet to Neumann b.c, was obtained previously by Kostov et al. in 
the context of two-dimensional gravity f48] . If ones adds strings (say L strings) , one has to choose 
if the leftmost string is blobbed or not. Clearly, if it is blobbed, it gives rise to half-loops when we 
go to open b.c (see figure [TH), so this case should be the same as the one of L — 1 strings when the 
leftmost one is unblobbed. This is in agreement with 

Kl{q)=Kl_, (82) 

in the case ri = '^^1^ ^ (see formula (|46l) ). If one wishes to write down a partition function on an 
annulus with (for example) open/ordinary b.c on each side, one should proceed as follows. First, one 
should note that to a unique configuration of the loops with open b.c correspond two configurations 
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with AS^ b.c, as shown in figure \T7\ One of these configurations contains a blobbed loop which 
winds around the annulus, the other does not. Then, to avoid double counting, we must count the 
blobbed loops which wind around the annulus with a fugacity 0. 




Figure 17: The two configurations for anisotropic special b.c with m — 1 leading to the open b.c 
configuration shown in figure [TBI (6). On the annulus (top and bottom of the lattice identified), one 
contains a blobbed loop (a) winding around the annulus, the other does not (b). 

This can be done easily, if we recall that the partition function has the Markov trace structure 
((53l) in the case with blobbed loops. The only thing one has to do is to modify the coefficients 
and 2?" in formula ([57)1 . One has "Dq = 1, "D* = 0, V' = Vi = n and the relation ((56)) leads to 

_ sm(i-L)7 sin(i+L)7 ^ rpj^- j-gjg^^ jqjj could also be obtained by specializing the results of 

sin 7 ' sin 7 j ir' o 

[Hll]. We end up with the partition function 

where ri — . The modular transform of ((83)) is computed once again by performing a Poisson 

resummation, which yields 

/2r^Y^ sin(7/2-^.2i(2p+i)^ 
^owo.. = ^ (84) 

The identification of the g-factor of the open b.c goven and the associated boundary entropy Sopen = 
— log gopen is straightforward 

(2\"' sin (7/2) 

90pen = - I . 1/2 - (85) 

(sm7sm(7/(7)) 

6.3 Open/open boundary conditions on the strip or the annulus 

The weight of the half-loops vi turns out to be unimportant because in the scaling limit the b.c 
renormaHzes to fixed b.c (ie vi = 1) as soon as vi > 0. In the case of two boundaries {eg on a strip 
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or on the annulus), however there is a non-trivial scaling behaviour of the model. Let us consider 
the partition function 



Z = Tr{ W 

boundary 1 



1 + S,.B^ 



n 

boundary 2 



1 + S,.B^^ 



bulk <ij> 



(86) 



where Bi and B2 are two boundary fields. Expanding this expression to get the loop expansion, 
one gets a model with half-loops attached on each boundary. A half-loop attached on the boundary 
1 gets a weight i^i, a half- loop attached on the boundary 2 gets a weight V2, and a half-loop with 
one end on each boundary has a weight 1^12. These parameters are related to the fields i?i and B2 
by 



1^1 



Si 



V2 



V\2 = B1.B2 



(87) 



Again, on each boundary the b.c should fiow towards fixed b.c as soon as Bi > and B2 > 0. 
The scaling limit is thus independent of the field strengths. However, it should still depend on the 
angle between the two fields n ^],'u^ 1 = rr^ - We can then proceed as in the case with only 

^ IIS2II ^ 

one boundary, by completing the lattice to get a configuration of the 0{n) model with ASu/ASo 
b.c with ni — n2 = I. The relation between these boundary conditions and the open/open case is 
exact only if — 1^2 — ^- In that case, a pair of half- loops attached on both boundaries on the 



open/open lattice gives a doubly blobbed loop with weight ni2 



When J^i 7^ 1 or 1^2 7^ 1 there 



is no exact mapping in the discrete setting between the configurations. However, because the only 
non-trivial parameter remaining in the scaling limit is the angle between the fields in the open/open 
case and the loop weight ni2 in the AS^/ASo case, the two models have to be equivalent in the 
scaling limit if 



ni2 



^12 
iyih'2' 



(88) 



In particular, this has the direct consequence that in the sector without strings, the Virasoro 
characteilill for open/ open b.c is (|40l) where the parameters ri2 and 1^12 are related by 



^Vi_V2 



+ 2 = 2 cos 



(89) 



which is a consequence of l(34|) with t\ 



r2 



7r/7 — 1 



The relation ([89]) was known by Cardy 07], 



who derived it by the following Coulomb gas argument! I. Recall that in the bulk the oriented loops 



get weights 



= ±«7 



depending on their orientation. On the open lattice, one counts the boundary 



monomers with weights a, 5, fi (see figure fT8|) and one parameterizes 



a — Pie 
5 = p2e' 



(3 



Pie" 
P2e~ 



(90) 



^^Note that in the calculation which led to Ko{q) (see section UTTIl the number of half-loops with weight U12 is 
always even, because they correspond pairwise to a loop with weight ni2. This calculation has to be modified if one 
is also interested in the case of odd number of half-loops. 

^^We thank J. Cardy for the permission to include this independent argument [4 7) . which generalizes his results 
for percolation [35| . 
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then the unoriented loops 


get the 


weig 


hts 1^1. 






a/3 






^2 = 


6fi 






Vl2 = 


aS 




Noting that 2 cos (7/2) = 


Vn + 2. 


, one 


has 2 



-'■7/2 



=.-^7/2 



2p? COS (7/2) 
2p2 COS (7/2) 



(91a) 

(91b) 
(91c) 



==Vn + 2, which is the same as relation 



o 

m 




Figure 18: Coulomb gas for the open/open b.c [47] (which generalizes the arguments of [35]). As 
usual, the loops are oriented: in the bulk they get a weight e**''' depending on their orientation. 
The boundary monomers get weight a, f}, S, fi, depending on the orientation. 

((89|) with 29 = Now, as in section 14.1.31 it is argued that pi and p2 do not contribute to the 

universal part of the boundary free energy, since these boundary conditions are expected to flow 
towards fixed b.c. Then the universal behaviour only depends on the phase factors e**^^ for each 
loop wrapping around the annulus, and the rest of the argument goes exactly as the calculation of 
the character Kq (formula l(40|) ) performed in section [4. 1.31 

Strings could be added using the results for the AS^/ASo case. One could also compute full 
partition functions instead of simple Virasoro characters, using the Markov trace decomposition 
([53]) . and treating the blobbed loops which wind around the annulus with care, following what we 
did in section [631 We do not want to go through this here, and we turn now to a simple application 
of the open/open b.c to the Ising model. 

6.4 Applications to crossing probabilities of Ising clusters on an annulus 

As a particularly simple application of our open/open b.c results, we compute some quantities 
related to the crossing probabilities of Ising spin clusters on the strip and the annulus. So far in 
this paper, the 0{n) model was formulated on the honeycomb lattice. When 71 = 1 we thus get the 
Ising model on the honeycomb lattice. This Ising model is dual to the Ising model on the triangular 
lattice, and this is the main point we use in what follows. The Ising spin cluster on the triangular 
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K > K, 




Figure 19: Ising model on an annulus (periodic boundary conditions on the left and right sides, 
free on top and bottom) as a function of the coupling K. We are interested in the cluster interfaces 
which cross the annulus from the bottom to the top. In these examples, there is no such interface 
for K > Kc, 2 interfaces for K = Kc and 6 for K < K^.. 



lattice have boundaries which are exactly the 0(1) loops above with a weight x per monomer. Note 
that X is related to the coupHng energy between two neighbouring spins of the triangular lattice 
E = -K(Ji(jj hy x = e^^. 

Let us consider the Ising model on an annulus of aspect ratio t — T/L where T is the periodic 
direction. We choose free boundary conditions for the spins. We are interested in the probability 
P{t) that there is at least one cluster boundary crossing the annulus from one boundary to the 
other (figure [H]) . Of course, such a probability behaves non-trivially with the coupHng K. When 
K > Kc, the model renormalizes towards a trivial theory with all spins frozen in a common state. 
Thus we expect P{t) = in that case. When K < Kc it renormalizes to a model where all the 
spins are independent, and they can be in either state with probability 1/2. This is nothing but 
critical percolation. The crossing probability for percolation cluster boundaries was first computed 
by Cardy in [35]. It can be written as a ratio of Dedekind eta functions 

ppcrc. ^ _ v{itMit/6)^ 

^ ^ ^ 77(zr/2)2,y(zr/3) ■ ^ ' 
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The case we are interested in is of course the case of the critical Ising model K — Kc, which should 
lead to a non-trivial probability PJ*'™s(t). To compute this we need the partition function of the 
critical Ising model on an annulus with free/free (ordinary /ordinary in the terminology we used 
above) b.c. We could of course compute this with the above formalism, by plugging 7 = -I, g = | 
into relation ([54)1 (recall also ^45^). This would yield the standard result [37tl38] 

^frco/froc(T) = + XiAq) (93) 

where Xr,s is the usual Rocha-Caridi character and q — e^^'^ as above. We could also have computed 
this partition function within the open/open b.c framework, because free b.c for the spins on the 
triangular lattice are equivalent open b.c for the cluster boundaries. We then have Zopen/Open — 
Zfi-oe/frec- The Markov trace structure of Zopen/Open leads to 

Z Open/ Open — Ko{q, ni2 = 1) + tcrms independent of ni2. (94) 

In this expression, it is very easy to subtract the configurations which do not contribute to the 
crossing probability: it is the same expression with ni2 — instead of ni2 ~ 1. Since the parameter 
ni2 only appears via ri2 in the character KQ{q, 1112), one ends up with 

Pl'-^r)^ '^ \ > (95) 

where we used (|40)) with ri2 — 1 (ni2 — 1) and ri2 = 3 (ni2 — 0). The relation ([95)) can be 
reformulated in terms of Dedekind eta function, to get it in a form very close to ((92)) 



i,i„ ^ r,iiT)r)iiT/ 12)^ 

" ^ ' rj{iT/2yr]{iT/6)' ^ ' 



Many other quantities could be computed. Let us mention but one of them before we turn to 
the numerical check of these results. Instead of computing the probability that there is at least one 
cluster boundary which crosses the boundary, one could be more precise and ask how many of them 
there are. One could derive the probability that there are fc > 1 clusters which cross the annulus 
as a function of the aspect ratio. For critical percolation this was done by Cardy in |35], and we 
could do it for the critical Ising model. This would turn out to involve more complicated formulae, 
but we can focus on a simpler quantity. Consider the limit when r ^ 1, then the annulus looks 
like a very long strip. We can then ask how many cluster boundaries are crossing the strip per unit 
length (on average). Again, for percolation this is a result of Cardy [49] : there are (on average) ^ 
cluster boundaries per unit length which cross the strip if we assume that the width of the strip is 
1. Let us derive the similar result for the critical Ising model here. 



In the limit t ^ 1 the partition function Zopen/Open behaves as 

Zopen/Openi'^12) - g'''-i2.^i2"-/24 (97) 
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where vi2 and ri2 are related by ([89)) with vi — V2 — ^- The average number of cluster boundaries 
per unit length is in this limit 

bound./ n 1 ^ / \ 
= — Ov 12 ^Og/j Open/ Openyl' 12) 

dvi2 
2g sin(7/2) 

which yields for the critical Ising model. 
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Figure 20: Probability that there is at least one contractible spin cluster crossing the annulus for 
different sizes. The simulations are made on a triangular lattice. The aspect ratio r is kept fixed and 
is equal to r = 2/^/3. The horizontal lines correspond to P = 0, P = Pcrit. ising and P = Pcrit. perco. 
from bottom to top. The vertical one is = Kc- 

Numerical check: We checked these results for the Ising model by Monte Carlo simulations. Note 
that it is very important to be right at the critical point when we compute crossing probabilities 
numerically, otherwise we are quickly attracted by another critical behaviour (critical percolation 
if K < Kc, frozen spins if K > Kc). Because simulations are done in finite size, there is a small 
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shift in the effective critical coupling around Kc 05- It is then a bit difficult to catch the right 
critical coupling, because of these finite size effects. To define the quantity we actually measure 
numerically and which is to be compared to our prediction for Pcrit. ising, we proceed as follows. For 
each aspect ratio r, we compute the probability P(t, K, N) for several K around Kc and different 
systems of size A'' x {tN). We plot P{t, K, N) as a function of the coupHng K for the different sizes 
TV X (tN), and find that the successive curves all intersect in a very small region at some effective 
critical coupling Kc{t,N). This procedure is shown in figure Eifor r = 2/V3. The value of the 
probability at this point is the quantity we compare with our analytic results. For K < Kc we 
also get a measure for the percolation clusters crossing probability (figure [20]). We plot the results 
obtained in this way in figure [2TJ For the average number of cluster boundaries crossing a strip per 
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Figure 21: Probability that there exists at least one contractible spin cluster crossing the annulus 
as a function of its aspect ratio. There are two possible non-trivial scaling limits. One is of course 
when K = Kc, then the crossing probability is given by ([96]) . The second one is the low-coupling 
limit K <C Kc, when the spin clusters look like critical percolation clusters (|92p . 

unit length, we simulated a system of aspect ratio r = 20. We find that there are 0.432 ± 0.002 
cluster boundaries crossing the annulus per unit length (to be compared with the analytical result 
V3/4 = 0.4330...). 

7 Conclusion 

In this article we have introduced a dilute version of the conformal boundary loop model, the study 
of which was initiated in [4J and pursued in [H [6]. For clarity, we now recall the main results 
obtained in this rather long paper. 
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7.1 Summary 

The dilute model is defined by ^ in terms of 0{n) type vector spins, and by ifTTj) in terms of 
an ensemble of self-avoiding loops. In the spin representation, its most important physical feature 
is the possibility of attracting preferentially a subset of ni indices (or, by duality, the remaining 
n — ni indices) towards the boundary. Correspondingly, in the loop representation, boundary 
touching loops come in two versions (blobbed and unblobbed) with respective weights n, = ni and 
n, = n — rii. 

An essential ingredient is the identification of two sets of integrable weights in a corresponding 
lattice model. These read l(24ll - l(25l) for the square lattice, and (fT2l) in the honeycomb limit. It 
should be noticed that the integrable weights impose particular weights for the boundary monomers 
of blobbed and unblobbed loops. 

The ordinary and special transitions (denoted Ord and Sp) in the 0{n) model are well studied 
in the literature, and correspond in our setting to the special case where blobbed and unblobbed 
loops are indistinguishable. The integrable points found here complete the picture by defining a 
pair of anisotropic transitions AS^ and AS,. The physical interpretation of AS^ is that blobbed 
loops are critically attracted towards the boundary (i.e., they stand at a special transition), whereas 
the unblobbed loops are repelled from the boundary (i.e., they stand at an ordinary transition). 
This interpretation is validated by the results of section l43l on the fractal dimensions of the contact 
sets of each loop type. Note that AS, is obtained from AS, by exchanging rii and n — rii. The 
complete phase diagram (in the physical region < ni < n) is shown in figure HI 

Using arguments of boundary conformal field theory (BCFT) and Coulomb gas we have identified 
the boundary-condition-changing (B.C.C) operators corresponding to the above transitions. In 
particular, we have ^ri,ri+i for {AS,\Ord), and ^n.n for (AS^lOrd). We have also obtained the 
corresponding Virasoro characters l(46|) in the sector with L non-contractible loops. There can 
be combined into conformally invariant partition functions (|57l) which encode the full operator 
spectrum of the theory and the corresponding multiplicities. 

The boundary entropies of the conformal boundary states identified in this work are given 
in (|58|) -([6T 1) . They are consistent with the g-theorem [42J and the proposed phase diagram (see 
figure H]). We have also identified the operators perturbing the isotropic points along the anisotropic 
direction, which are <I>3,i at Ord and $3,3 at Sp. This shows that the anisotropic transition lines 
of figure |4] form a cusp where they join in Sp. Starting at the anisotropic point ^5*,, the operator 
perturbing in the unstable direction was identified in section [5741 as $1,3. 

Finally, we have shown in section [6] that open boundary conditions for the 0{n) model are a 
special case (ni = 1) of those considered here. They correspond to imposing a magnetic field on the 
boundary spins. As a simple geometrical application, we have derived a new crossing probability 
l(96|) for Ising domain walls, and found the number of cluster boundaries per unit length that cross 
an infinite strip. 

7.2 Outlook 

Following [6] one may consider a more general model on the annulus with two distinguished bound- 
aries and loop weights (|30l) - l(3T|) and fM)) . In this model, one can also introduce separate weights 
for contractible and non-contractible loops. We have here only alluded briefiy to this generalization, 
which will be treated more fully elsewhere [33] . 

While this work completes the conformal part of the program set out in [4], we have omitted 
here a whole range of algebraic questions. In particular, one may study the (rational) restrictions of 
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the conformal boundary loop models when some of the parameters take particular "magic" values. 
This will also be treated in [33j . 

Another future direction would be to exploit the integrable solutions (|24l) - l(25|) of the reflection 
equation to set up the corresponding Bethe ansatz equations. Presumably this would put our results 
for the critical exponents and the spectrum generating (partition) functions on a rigorous basis, 
and would allow to deduce non-universal quantities such as the surface free energies. 

Also, let us discuss what may happen when boundary loops can be decorated with k different 
orthogonal blobs. The weight of an ?-blobbed loop (for i — l,2,...,fc) is and its boundary 
monomers have fugacity Wi. Obviously, n = J2i=i '^ji ^^d the preceding discussion corresponds to 
the case k = 2. When several Wi coincide, the corresponding blobbed loops are indistinguishable 
and their weights rii may be regrouped. In particular, when all Wi coincide we recover the transitions 
Ord and Sp. 

Suppose now that the set {wi} takes precisely two different values, with fci weights equal to w 
and k — ki weights equal to Hj, and w > w. Regrouping the corresponding loop weights we are 
then, in fact, in the situation k = 2 where the previous results apply. Correspondingly, we have 
an anisotropic special transition, with ki groups of indices standing at a special transition, and the 
remaining k — ki groups standing at an ordinary transition. In the full parameter space {wi}, this 
critical point has ki unstable directions (corresponding to moving a Wi away from w) and k — ki 
stable directions (corresponding to moving a Wi away from iB). This argument yields a total of 2*^ 
critical points (including Ord and Sp). We cannot exclude the existence of further multi-critical 
points with finite values of {wi}, but we conjecture that no such point exists. 

Finally, we expect that there are applications of this work to other physical problems. One of 
these is the re-intepretation of our conformal boundary conditions in terms of boundary degrees 
of freedom critically coupled to the bulk. While such a reinterpretation involves Uq[su{2)] spins in 
the case ri an integer, it is not clear what is the meaning of ri not integer, and whether this has 
something to do with sZ(2,IR). One could also wonder whether our CBCs have any relation with 
boundary bound states in the S matrix description of the bulk CFTs [50]. Finally, we note that the 
Hmit n ^ 2 is deeply related with the Kondo model. What happens to our phase diagram in this 
case is also an open problem. 
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